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Abstract 

An  expression  for  the  semi-major  axis  as  a  function  of  time  may  be  determined  by 
performing  a  series  reversion  and  inversion  of  Lambert’s  Time  Function.  Since  the 
resulting  series  contains  a  singularity,  it  is  desirable  to  perform  only  a  reversion  on  the 
original  series  to  obtain  an  expression  for  the  inverse  of  the  semi-major  axis.  Using  a 
Lagrange  expansion  to  obtain  the  coefficients  for  this  series  is  very  computer  intensive. 
Therefore,  alternative  methods  are  presented.  Also,  each  series  was  expanded  into  a 
continued  fraction  which  provided  greater  accuracy  than  the  series  using  the  same  num¬ 
ber  of  coefficients.  The  accuracy  was  found  to  be  dependent  upon  the  number  of  series 
coefficients  used,  the  transfer  time,  and  the  transfer  angle. 
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A  NUMERICAL  INVESTIGATION  OF  THE  SERIES  REVERSION/INVERSION  AND 
SERIES  REVERSION  OF  LAMBERTS  TIME  FUNCTION 


I.  Introduction 

One  of  the  classic  problems  in  orbital  mechanics  is  the  determination  of  orbital 
parameters  using  two  position  vectors  and  the  time  of  flight  of  a  satellite.  In  1761, 
Lambert  developed  a  theorem  (later  proved  by  Lagrange)  which  states  the  time  to  tra¬ 
verse  an  elliptic  arc  (also  true  for  any  conic  trajectory  (2:/  ))  is  dependent  only  on  the 
semi-major  axis,  the  perimeter  of  the  triangle  formed  by  the  initial  and  final  position 
vectors,  and  the  chord  joining  the  initial  and  final  points  of  the  arc.  The  geometry  is 
shown  in  Figure  1.  Throughout  this  research  the  radius  vectors  were  held  constant, 
r,  -  r2-  1  au  and  the  transfer  angle  was  varied  from  10  to  360  degrees. 


The  significance  of  Lambert’s  theorem  can  be  seen  in  Figure  2.  Consider  the 
elliptic  arc  joining  the  fixed  points  P  and  Q.  According  to  Lambert’s  theorem  the  shape 
of  the  transfer  ellipse  from  P  to  Q  may  be  altered  without  changing  the  time  of  flight 
by  moving  the  foci  F  and  F"  and  keeping  the  semi-major  axis  constant.  Therefore,  the 
focus  F  may  be  moved  to  F,  and  the  focus  F'to  F\  without  changing  the  time  of  flight. 
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As  the  focus  F  is  moved  counterclockwise  and  the  focus  p'\s  moved  clockwise,  the 
ellipse  becomes  very  flat  until  the  foci  are  at  F  2  and  F  \  .  The  time  of  flight  has  not 
changed  and  the  orbit  is  now  rectilinear.  Thus,  the  time  /  may  be  computed  by  elemen¬ 
tary  methods. 


Figure  2.  Geometrical  Significance  of  Lambert’s  Theorem  (2:72) 


To  obtain  Lambert’s  Time  Function  the  energy  equation  for  a  two-body  orbit  may 
be  written  as 


(1) 


where  r  is  the  radial  distance  between  the  two  bodies,  p  is  the  gravitational  constant,  a  is 
the  semi-major  axis,  and  t  is  the  time.  Since 


QF2  +  PF2  -  r,  +  rz  and  QF2-PF2-c 


then 


PF2ms-c  and  QF2-s 


where  s  -  (r ,  +  r 2  +  c)/2  .  Integrating  Eq  (1)  produces 

1  I"  rdr 
■J\i^s-c^2r-rz/a 

Employing  the  following  change  of  variable: 


(2) 


2 


r  -  a(  1  -  cos$) 


Eq  (2)  may  be  written  as 


where 


—  ( 1  -  cos4>)c/$ 

P 


a  =  2sin 


and 


3 


This  produces 


t 


[(o-sina)-((5-sinP)] 


(3) 


This  is  Lambert’s  Time  Function  for  an  elliptic  trajectory  with  a  transfer  angle  less  than 
n  and  a  flight  time  less  than  the  minimum  energy  transfer  time.  A  minimum  energy 
transfer  is  an  elliptical  transfer  where  a  =  s/2  .  It  may  also  be  shown  for  a  transfer 
angle  greater  than  rt  Lambert’s  Time  Function  is  given  by 


( 


— [(a -sin  a)  +  (3 -sin  (i)] 


(4) 


This  form  of  Lambert’s  Time  Function  is  difficult  to  use  since  the  transfer  time  is 
often  known  and  the  semi-major  axis  is  the  parameter  that  needs  to  be  determined. 
Since  the  equation  is  transcendental  in  the  semi-major  axis  root  finding  techniques  are 
employed  in  order  to  determine  the  semi-major  axis.  A  method  of  expressing  the  semi¬ 
major  axis  as  a  function  of  the  transfer  time  is  desired. 


With  the  use  of  hypergeometric  functions,  Lambert’s  Time  Function  for  an  ellipti¬ 
cal  trajectory  may  be  written  as  an  infinite  series: 


where 


3 


is  the  parabolic  transfer  time.  This  is  the  flight  time  of  an  object  traveling  on  a  para¬ 
bolic  trajectory  between  two  position  vectors  where  the  origin  of  the  position  vectors  is 
the  focus  of  the  parabola.  Lambert’s  Time  Function  for  a  hyperbolic  trajectory  may  also 
be  expressed  as  an  infinite  series  identical  to  Eq  (5)  except  in  the  sign  of  the  argument, 
i.  e.,  (-^)  ,  and  the  series  coefficients  are  identical. 

A  series  reversion  and  inversion  has  been  performed  on  Eq  (5)  in  order  to  obtain 
the  semi- major  axis  as  a  function  of  time  (4:9-13): 

(  ^  j  -  S,7"‘  +  B2T°  +  B3Tl  +  S472  +  ...  (6) 

The  coefficients  for  this  series,  referred  to  as  the  RI  series,  are  obtained  using  a  numer¬ 
ically  efficient  matrix  formulation.  One  purpose  of  this  research  is  to  investigate  the 
accuracy  of  the  RI  series.  Due  to  the  existence  of  a  singularity  in  the  RI  series  at  T  -  0 
it  is  desirable  to  perform  only  a  series  reversion  on  Eq  (5)  to  obtain  the  inverse  of  the 
semi-major  axis  as  a  function  of  time: 

^yaj+azT^a^+aj^..,  (7) 

This  series  will  be  referred  to  as  the  R  series.  Several  methods  of  obtaining  the  coeffi¬ 
cients  for  the  R  series  are  presented  as  well  as  the  accuracy  of  each  method. 
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II.  Analytical  Development 


Measuring  the  Accuracy 

The  accuracy  of  each  method  was  evaluated  in  the  following  manner:  1)  For  a 
given  geometry  and  transfer  time  a  value  for  the  semi-major  axis  was  obtained  using  Eq 
(6)  or  (7)  for  a  given  number  of  terms;  2)  The  time  was  calculated  from  Eq  (3)  or  (4) 
using  the  value  for  semi-major  axis  obtained  in  step  1;  3)  The  absolute  value  of  the 
difference  between  the  calculated  time  and  the  true  time  was  then  calculated.  The  value 
obtained  in  step  3  is  called  the  Time  Residual  and  was  employed  as  a  measure  of  the 
accuracy  of  each  method. 

Calculating  the  RI  Series  Coefficients 

To  calculate  the  coefficients  for  the  RI  series,  multiply  Eq  (6)  by  to  obtain: 

T-s'(sM^h  (8> 

Substituting  the  series  solution  for  T  given  by  Eq  (5)  into  Eq  (8)  yields: 


Now  differentiate  Eq  (9)  with  respect  to  (^)  to  obtain: 


Evaluating  Eq  (10)  at  (£)  -  0  yields: 

4, -S, 

Taking  three  more  derivatives  and  evaluating  each  at  (^)  -  0  produces  the  following: 
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(11) 


A2  -  At32 
a3  -  a2bz*  a\b2 

A 4  ”  A 3 B 2  +  2  A  j  /1 2 3  +  /t  ^ B 4 

Each  successive  derivative  becomes  increasingly  more  complex.  Solving  Equation  set 
(11)  for  the  B,  coefficients  yields: 


di 

d, 

axa3-a\ 

a! 

A*  A  4  -  3  A  {  A 2  A  3  +  2  A% 


Equation  set  (12)  may  be  expressed  in  matrix  form  as: 


(12) 


(13) 


A  method  of  generating  the  Q  matrix  in  Eq  (13)  has  been  discovered  (4:12). 
Knowing  the  elements  of  the  first  two  rows  of  the  matrix  all  of  the  remaining  elements 
may  be  obtained  using  the  following  procedure: 

.  i  Ai  i 

1-  Qu  =  t,  ■  q2'°  ~  7]  '  =  7, 


2.  q„  =  (Q,-i  i  .  q.-z  i . q  n  ) '  (q  i 1  .  qfz. 1 . q,- i,-i).  i*j>2 

3.  q„  =  -  j-<(.q,2.q,i . qu)-(.A2.A3 . a,)} 

n  1 

4.  S,.|  =  (q„  .  q,2 . q„)-(A2.  A3 . /!,.,) 
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The  first  element  of  each  row  must  be  calculated  last,  since  it  is  dependent  upon  all  of 
the  other  elements  in  that  row.  For  example,  to  calculate  fi4  first  form  element  q32 : 

Q33  =(g2i  ■Qu)  *  (Q  12  •  c/zz) 


(- 


A3> 


then  q 32  : 


Q 32  _(<72ii£7ii)‘(<7ii'<72i) 


then  q3,  : 


Q  31  ~  ~  A  {(<732'<733)’('^2,/^3)) 

n  I 

1  /  2A\  /t3\  2A22-AiA3 

=  ~  A,  \  /It  +  A3J  A ? 


Now  that  all  of  the  elements  of  the  third  row  are  known,  fi„may  be  calculated  by  taking 
the  dot  product  shown  in  step  4.  Each  successive  coefficient  may  be  calculated  by 
repeating  steps  2,  3,  and  4.  The  first  26  coefficients  for  a  transfer  angle  of  60  degrees 
are  listed  in  Table  1. 
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TABLE  1:  RI  Series  Coefficients  for  0-60 


n 

Bn 

1 

3.4766271 0943897  IE-001 

O 

4. 

5.601954800296930E-001 

3 

1 .569332354004 1 66E-001 

4 

-2.072 127624535575E-002 

5 

1 .34 1 48508098 1 992E-002 

6 

-9.59 14691 7620452  IE-003 

7 

5.932190937546468E-003 

8 

-3. 15903 1868854800E-003 

9 

1 .491 6643 14059560E-003 

10 

-6.607 1 178240781 35E-004 

11 

3.008808653647832E-004 

12 

-1.524580592486835E-004 

13 

8.421050704976452E-005 

14 

-4.62684562 1302778E-005 

15 

2.360695540204460E-005 

16 

- 1.1031  1 34086 18 128E-005 

17 

4.9447807555225 15E-006 

18 

-2.4088 10757970059E-006 

19 

1.33533831 1037049E-006 

20 

- 1 .4282095435 191 88E-006 

21 

-2. 1608448 15080054E-006 

22 

2. 1405 1 36067 1 874 1 E-005 

23 

- 1. 1608 19447250 105E-004 

24 

6.927 153735887259E-004 

25 

-3.255901 276 133955E-003 

26 

1. 38257488 1426990E-002 

Calculating  the  R  Series  Coefficients 

The  coefficients  for  the  R  series  may  be  obtained  by  using  the  Lagrange  expansion 
theorem  (1).  Letting  x-(^)  ,  Eq  (5)  may  be  written  in  the  form: 

x-x0  +FKx)  (14) 


with  v0  -  0  ,  where 


♦  (*) 


_ 1 _ 

4,  +  A2x*  A3x2*... 


(15) 


and  the  An  terms  are  the  coefficients  in  Eq  (5).  Now  x  may  be  expressed  as  a  power 
series  in  T  using  the  Lagrange  expansion  theorem  with  the  coefficients  given  by 


a 


n 


1  d"' 
n  1  dx"'1 


♦  (*)" 


X  *  X 


0 


8 


The  following  recursive  algorithm  was  used  to  calculate  the  coefficients  of  Eq  (7) 
(3:212-215).  Define 

d k 

D'k  =  — -  4>  ( -V ) "  for  A:  =  0 ,  1 . n-1 

dx 


and 


dxk 


k<*> 


(*„)-♦ 


(*> 


The  algorithm  now  becomes 


,<*>_  1  y  k>  , 


(*- O 


V'f  1  ^ 


for  *-1.2. 

beginning  with  the  initial  values 


D'o-W-- 


1_ 


The  coefficients  can  now  be  calculated: 


a  -  —  D"  i 
•  n'  "  1 


Algebraically,  the  first  five  coefficients  are 


. .  n  -  1 


1 

'*  "  A, 


A  3  A\ 

- -  +  2  — 

A\  ** 


Ak  A  2  A  3  d  2 

- :  +  5 I 5  — 

A*  A t 


/l5  6/t2/l4  +  3/l§_2i  /1 2^3  ,  ,  A\ 


A* 


A7, 


A' 


■*  14-J 


(16) 
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Each  successive  coefficient  becomes  increasingly  more  complex.  Table  2  lists  the  R 


series  coefficients  for  a  transfer  angle  of  60  degrees. 

TABLE  2:  R  Series  Coefficients  for  8-60° 


n 

a* 

1 

2.876351039445733 

2 

-4.634718652746828 

3 

6.169638843024049 

4 

-7.677730891005001 

5 

9.199104623157812 

6 

-10.731091604324702 

7 

12.266169580346382 

8 

-13.800785808657999 

9 

15.334487956965681 

10 

-16.867784047477109 

1 1 

18.401065122700818 

12 

-19.934436495815678 

13 

21.4678641 1 1717599 

14 

-23.001301593856162 

15 

24.534730221361581 

16 

-26.068150929596808 

17 

27.601569133880563 

18 

-29.134987387096338 

19 

30.668404962779338 

20 

-32.201822741326474 

21 

33.735279494304024 

22 

-35.269089203411149 

23 

36.804909127457499 

24 

-38.350692931914622 

25 

39.939922586593518 

26 

-41.694275534386783 

Since  the  recursive  algorithm  used  to  calculate  the  R  series  coefficients  is  slow  in 
terms  of  computer  time  a  more  efficient  method  is  desirable.  The  algebraic  form  of  the 
coefficients  were  formulated  into  a  matrix  in  order  to  determine  whether  a  pattern 
existed  such  as  for  the  RI  series.  Alternate  methods  of  obtaining  the  R  series 
coefficients  examined  the  behavior  of  the  coefficients  as  functions  of  n  and  0. 
Additionally,  each  series  was  expressed  as  a  continued  fraction  so  that  fewer  coefficients 
would  be  required  in  order  to  achieve  the  same  accuracy  as  the  series. 
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Matrix  Method 


Equation  set  (16)  may  be  expressed  in  the  following  matrix  form: 


Since  the  first  row  and  the  first  column  of  P  contain  all  zero  elements,  except  the  first 
element,  we  can  look  at  the  submatrix  P 


P 


_1_ 

A\ 


_1_ 

Ax 

2A2 
’  A* 

5  A\ 

~aT 

14/^2 

t- 


_i_ 

At 

5A2 

~^T 

21  A\ -3/1,  A  3 


0  0 

0  0 

*  ° 

6A2  J_ 

Al  A ? 


0 

0 


0 

0 


R 


0 

1 


A\ 


i  J 


The  first  column  of  the  R  matrix  has  the  general  form: 

(-  1),-|(2Q'  A‘2' 
""  i»(i+l)'  Af' 


and  the  diagonal  has  the  general  form: 


II 


For  i  >2,  the  element  to  the  left  of  the  diagonal  element  has  the  general  form: 


R 


t.  ('  t 


0  +  2  )A2 
A\" 


These  were  the  only  general  terms  that  were  ascertained.  Further  investigation  may 
reveal  a  pattern  that  would  be  beneficial  in  order  to  generate  the  entire  matrix. 


Approximating  the  R  Series  Coefficients 

Two  methods  of  approximating  the  R  series  coefficients  were  investigated.  The 
first  method  (Method  A)  involved  obtaining  the  coefficients  as  a  function  of  n.  The 
coefficients  were  plotted  versus  n  for  each  transfer  angle.  A  least  squares  analysis  was 
then  employed  to  fit  a  polynomial  expressing  a„as  a  function  of  n  to  each  curve.  The 
second  method  (Method  B)  involved  obtaining  the  coefficients  for  the  R  series  as  a 
function  of  the  transfer  angle.  Each  coefficient,  a,  through  a,6,  was  plotted  versus  the 
transfer  angle.  Once  the  plots  were  obtained,  a  least  squares  analysis  was  performed  to 
obtain  a  polynomial  in  the  transfer  angle  for  each  a„. 


Continued  Fraction  Expansion 

The  final  method  examined  involved  expanding  the  series  into  a  continued  frac¬ 
tion.  This  method  was  applied  to  both  the  R  series  and  the  RI  series.  A  continued 
fraction  typically  provides  the  same  accuracy  as  the  series  but  uses  fewer  coefficients. 
The  general  form  of  a  continued  fraction  is 


In  general,  the  function 

,,  „  C|0  +  CUX  +  CI2X2+  .  .  . 

fix) - - - 

Coo  +  C0|X+CozX‘:+  .  .  . 

may  be  expanded  into  the  continued  fraction 


(17) 
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q  .  ^_10  C20  X  C  30  X 


The  coefficients  clk  are  computed  recursively  by  the  formula 


C 1-2.0  c i-2. k’  I 


',-1.0  L  I  .*•  1 


where  j  >  2  . 


The  RI  series  may  be  written  as 


and  the  R  series  may  be  written  as 


2  a  B , 
T-7+/(r) 


where 


/(T)_c10  +  cmT  +  c12T2  +  . 


For  the  RI  series  c,»  =  fi*.2  and  for  the  R  series  cu-at>1  for  A: -0.1,2 . 

Eq  (22)  is  of  the  form  of  Eq  (17)  where 

c00  -  1  and  cQk  -  0  for  Ac-1, 2.  3... 

The  coefficients  for  the  continued  fraction  expansion  may  now  be  calculated  using  Eq 
(19).  If  the  continued  fraction  is  terminated  so  that 


B  k  c ,  n  c w 7  C30T 

f(T)  _  _i  _  0;  -12  ,_?2_,_2£_ 

Qk  L  c00  CI0  C20 


then  ^  is  termed  the  Acth  convergent  of  the  continued  fraction.  The  convergents  are 
found  using  the  following  formula; 

B k  _  <*kBk- 1  +  bkPk-z 
Qk  a*Q«-i  +  bkQk.2 

where 


'■'*-1.0 
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“  cto^ 


k>  1 


and 


b 


c 


10  ' 


P0m  0.  P|-6[.  Q0“l.  and  Q,-l 

Algebraically,  the  first  four  convergents  are: 

P  i  f*2  _  cfo 

Qi  c10-cnr 


f’o  c  10C  I  I  (C  IqC  12  cf ,  )7~ 

Qa  ~  c  n  +  c  12  T 

f  4  _  (c  I0C  1 1  _cmc12)+(_2C|0CliC12  +  Cn  +  c]0C  |3)7 

Q4  (cN  -c10c12)  +  (cl0c13-c11c,2)7'  +  (cf2-cl,C|3)72 


(25) 


Equation  set  (25)  shows  the  &th  convergent  requires  up  to  the  coefficient  c,.*.,  . 

This  means  for  the  RI  series  the  coefficients  up  to  Bk. ,  are  required  and  for  the  R  series 
the  coefficients  up  to  a*  are  required.  The  accuracy  of  the  Arth  convergent  was  com¬ 
pared  to  the  accuracy  of  each  series  using  k  terms. 
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III.  Numerical  Investigation 


In  this  section  all  of  the  figures  shown  are  for  a  transfer  angle  of  60  degrees.  All 
other  transfer  angles  examined  were  similar  and  representative  samples  are  shown  in  the 
indicated  appendix.  To  measure  the  accuracy,  the  non-dimensional  time,  T  (defined  in 
Eq  (5)),  was  varied  from  -.95  to  1.00. 

Accuracy  of  RI  and  R  Series 

Using  the  method  described  in  the  previous  section,  the  accuracy  of  the  RI  and  the 
R  series  was  evaluated. 

Figures  3a  and  3b  show  plots  of  the  Time  Residual  versus  T  for  a  transfer  angle  of 
60  degrees  using  5  and  15  terms  of  the  RI  series.  Plots  for  other  transfer  angles  are  in 
Appendix  A.  Note  the  lack  of  a  data  point  at  T  -  0  .  This  is  due  to  the  singularity 
caused  by  the  first  term  of  the  series.  Five  terms  of  the  series  provides  good  accuracy 
for  j-T i  <  .50  while  15  terms  provides  good  accuracy  over  the  interval  |T|  <  .85.  The 
absolute  value  of  the  maximum  Time  Residual  for  transfer  angles  from  30  to  360 
degrees  is  plotted  in  Figures  4a  -  4c  for  various  intervals  of  time.  As  IT!  -»  1  the  maxi¬ 
mum  absolute  value  of  the  Time  Residual  increases  for  a  constant  number  of  terms. 
Therefore,  as  i  7*  I  -*  1  additional  terms  are  required  to  maintain  a  specified  degree  of 
accuracy.  For  |T|  <  .25  15  terms  provides  one  to  two  orders  of  magnitude  greater  accu¬ 
racy  than  5  terms.  For  \  T\<  .50  15  terms  provides  approximately  three  orders  of  magni¬ 
tude  greater  accuracy  than  5  terms.  For  |T|  <  .85  15  terms  provides  two  to  three  orders 
of  magnitude  greater  accuracy  than  5  terms  up  to  about  240  degrees.  Beyond  240 
degrees  there  is  only  about  one  order  of  magnitude  difference  between  5  and  15  terms. 
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TIME  RESIDUAL 


Rl  SERIES,  5  TERMS 
THETA  =  60  DEGREES 


1.0E-002  n 


4.2E-017  H 


-1.0E-002  : 

<  — 2.0E-002  : 

g  : 

£  — 3.0E-002  : 

I  -4.0E-002  : 

-5.0E-002  : 


— 6.0E-002  H 


*”7.0E— 002  t  i  t  "t  t  |  i  i  i  i — prr  r-r-]  n— t-tjttttti  r  r  t  y  r  t  t  t  t  T~~r~r— r~i 

-1.00  -0.75  -0.50  -0.25  0.00  0.25  0.50  0.75  1.00 

TIME  (t/tp  -  1) 


Figure  3a.  Time  Residual  vs.  T  for  a  Transfer  Angle  of  60°  Using  5 

Terms  of  the  RI  Series 


2.5E-005  n 
1.6E-019  : 
-2.5E-005  : 
-5.0E-005  : 
-7.5E— 005  : 


Rl  SERIES.  15  TERMS 
THETA  =  60  DEGREES 


-1.0E-004  : 

-1.3E— 004  : 

-1.5E-004  ^  • 

—  1.6E— 004  I  i  l  l  l — I  i  I  i  l  I  i  i — i  l  I  !  I  l  l  |  i — r— r- T— j— i— m — !  i  n— I  r'  “n — r— ] 

-1.00  -0.75  -0.50  -0.25  0.00  0.25  0.50  0.75  1.00 

TIME  (t/tp  -  1) 

Figure  3b.  Time  Residual  vs.  T  for  a  Transfer  Angle  of  60°  Using  15 

Terms  of  the  RI  Series 
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Rl  SERIES 


T  =  -.25  to  .25 


0  50  60  90  120  150  180  210  240  270  300  330  360 

ANGLE  (DEGREES) 


Figure  4a.  Maximum  Absolute  Value  of  the  Time  Residual  vs.  Trans 

fer  Angle  Using  the  RI  Series  for  T  -  -  .25  to  .25 


0  30  60  90  120  150  180  210  240  270  300  330  360 

ANGLE  (DEGREES) 


Figure  4b. 


Maximum  Absolute  Value  of  the  Time  Residual  vs.  Trans 
fer  Angle  Using  the  RI  Series  for  T  -  -  .50  to  .50 


0  30  60  90  120  150  180  210  240  270  300  330  360 

ANGLE  (DEGREES) 

Figure  4c.  Maximum  Absolute  Value  of  the  Time  Residual  vs.  Trans¬ 

fer  Angle  Using  the  RI  Series  for  7  -  -  .85  to  .85 

Plots  of  Time  Residual  versus  T  for  a  transfer  angle  of  60  degrees  using  5  and  15 
terms  of  the  R  series  are  shown  in  Figures  5a  and  5b.  Plots  for  other  transfer  angles  are 
in  Appendix  B.  The  accuracy  of  the  R  series  behaves  similarly  to  the  RI  series.  How¬ 
ever,  5  terms  provides  good  accuracy  for  |7|  <  .25  which  is  a  much  smaller  interval  than 
the  RI  series.  Fifteen  terms  also  provides  good  accuracy  over  a  smaller  interval  than  15 
terms  of  the  RI  series:  i  7"  j  <  .50.  Figures  6a  and  6b  show  plots  of  the  absolute  value  of 
the  maximum  Time  Residual  versus  transfer  angle  for  various  time  intervals.  For 
[7|  <  .25  15  terms  provides  about  four  orders  of  magnitude  greater  accuracy  than  5 
terms.  However,  the  R  series  is  not  as  accurate  as  the  RI  series  utilizing  the  same 
number  of  terms  and  has  a  smaller  interval  of  convergence. 


18 


ANGLE  (DEGREES) 

Figure  6a.  Maximum  Absolute  Value  of  the  Time  Residual  vs.  Trans 

fer  Angle  Using  the  R  Series  for  T  -  -  .25  to  .25 


ANGLE  (DEGREES) 

Figure  6b.  Maximum  Absolute  Value  of  the  Time  Residual  vs.  Trans 

fer  Angle  Using  the  R  Series  for  T  -  -  .50  to  .50 
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Method  A:  R  Series  Coefficients  as  a  Function  of  n 

Figure  7  shows  the  absolute  value  of  a„is  approximately  a  linear  function  of  n  for 

n-  1  to  approximately  n-  25.  A  least  squares  analysis  was  performed  using  20  terms 
and  the  first  26  terms  for  a  transfer  angle  of  60  degrees  in  order  to  obtain  a  polynomial 
equation.  The  residual  shown  in  Figures  8a  and  8b  is  the  difference  between  the  value 
of  a  „  calculated  using  the  polynomial  and  the  true  value.  In  Figure  8a  the  polynomials 
were  calculated  using  the  first  20  terms  and  in  Figure  8b  the  polynomials  were  calcu¬ 
lated  using  the  first  26  terms.  Using  20  terms  resulted  in  slightly  smaller  residuals.  In 
order  to  reduce  the  residuals  for  the  first  few  terms  two  polynomials  were  found.  The 
first  polynomial  was  for  n  -  1  to  n-m*  1  where  m  is  the  degree  of  the  polynomial. 
The  second  polynomial  was  for  n-m  +  2  to  n  -  20 . 


n 


Figure  7.  Absolute  Value  of  a„vs.  n  for  a  Transfer  Angle  of  60° 
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Figure  8a.  Absolute  Value  of  the  Coefficient  Residual  vs.  n  for  a 

Transfer  Angle  of  60°  Using  20  Terms 


Figure  8b.  Absolute  Value  of  the  Coefficient  Residual  vs.  n  for  a 

Transfer  Angle  of  60'  Using  26  Terms 


To  determine  which  degree  of  polynomial  would  provide  the  most  accurate  results, 
the  absolute  value  of  the  "a  Residual"  (defined  as  the  difference  between  the  semi-major 
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axis  calculated  using  the  true  coefficients  and  the  value  calculated  using  the  polynomial 
approximation)  was  plotted  versus  T.  In  order  to  plot  the  absolute  value  of  the  a 
Residual  on  a  logarithmic  scale  two  plots  were  required:  one  for  positive  values  of  T 
and  the  other  for  negative  values  of  T.  Polynomials  of  degree  two  through  degree  nine 
were  evaluated.  Figures  9-11  show  these  plots  using  5,  10,  and  15  terms  respectively. 


0.00  0.10  0.20  0.30  0.40  0.50  0.60  0.70  0.80  0.90  1.00 
TIME  (t/tp  -  1) 


Figure  9a.  Absolute  Value  of  the  a  Residual  vs.  Negative  T  for  a 

Transfer  Angle  of  60°  Using  5  Terms  with  the  Coefficients 
Calculated  Using  Polynomials  of  Degree  2-5 
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0.00  0.10  0.20  0.30  0.40  0.50  0.60  0.70  0.80  0.90  1.00 
TIME  (t/tp  -  1) 


Figure  9b.  Absolute  Value  of  the  a  Residual  vs.  Positive  T  for  a 

Transfer  Angle  of  60°  Using  5  Terms  with  the  Coefficients 
Calculated  Using  Polynomials  of  Degree  2-5 


0.00  0.10  0.20  0.30  0.40  0.50  0.60  0.70  0.80  0.90  1.00 
TIME  (t/tp  -  1) 


Figure  9c.  Absolute  Value  of  the  a  Residual  vs.  Negative  T  for  a 

Transfer  Angle  of  60°  Using  5  Terms  with  the  Coefficients 
Calculated  Using  Polynomials  of  Degree  6-9 
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0.00  0.10  0.20  0.30  0.40  0.50  0.60  0.70  0.80  0.90  1.00 
TIME  (t/tp  -  1) 


Figure  9d.  Absolute  Value  of  the  a  Residual  vs.  Positive  T  for  a 

Transfer  Angle  of  60s  Using  5  Terms  with  the  Coefficients 
Calculated  Using  Polynomials  of  Degree  6-9 
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Figure  10a.  Absolute  Value  of  the  a  Residual  vs.  Negative  T  for  a 

Transfer  Angle  of  60°  Using  10  Terms  with  the  Coeffi¬ 
cients  Calculated  Using  Polynomials  of  Degree  2-5 
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Figure  10b.  Absolute  Value  of  the  a  Residual  vs.  Positive  T  for  a 

Transfer  Angle  of  60°  Using  10  Terms  with  the  Coeffi¬ 
cients  Calculated  Using  Polynomials  of  Degree  2-5 
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Figure  10c.  Absolute  Value  of  the  a  Residual  vs.  Negative  T  for  a 

Transfer  Angle  of  60°  Using  10  Terms  with  the  Coeffi¬ 
cients  Calculated  Using  Polynomials  of  Degree  6-9 
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Figure  lOd.  Absolute  Value  of  the  a  Residual  vs.  Positive  T  for  a 

Transfer  Angle  of  60°  Using  10  Terms  with  the  Coeffi¬ 
cients  Calculated  Using  Polynomials  of  Degree  6-9 
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Figure  11a.  Absolute  Value  of  the  a  Residual  vs.  Negative  T  for  a 

Transfer  Angle  of  60' Using  15  Terms  with  the  Coeffi¬ 
cients  Calculated  Using  Polynomials  of  Degree  2-5 
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Figure  lib.  Absolute  Value  of  the  a  Residual  vs.  Positive  T  for  a 

Transfer  Angle  of  60°  Using  15  Terms  with  the  Coeffi 
cients  Calculated  Using  Polynomials  of  Degree  2-5 
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Figure  11c.  Absolute  Value  of  the  a  Residual  vs.  Negative  T  for  a 

Transfer  Angle  of  60°  Using  15  Terms  with  the  Coeffi 
cients  Calculated  Using  Polynomials  of  Degree  6-9 
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Figure  lid.  Absolute  Value  of  the  a  Residual  vs.  Positive  T  for  a 

Transfer  Angle  of  60°  Using  15  Terms  with  the  Coeffi¬ 
cients  Calculated  Using  Polynomials  of  Degree  6-9 


Figures  9-11  indicate  the  third  degree  polynomial  is  the  most  accurate  approxi¬ 
mation  over  the  time  interval  T  -  -  .95  to  1.00  .  Figures  10b  and  lOd  show  that  at 
approximately  T  -  .80  all  of  the  polynomials  give  poor  results  using  10  terms. 

Checking  the  a  Residual  using  6,  8,  12,  and  18  terms  showed  a  similar  spike  at  the  same 
value  of  T.  Therefore,  when  using  this  method  of  approximating  the  R  series  coeffi¬ 
cients  only  an  odd  number  of  terms  should  be  used. 

The  accuracy  of  this  approximation  using  the  third  order  polynomial  was  analyzed 
in  the  same  manner  as  the  true  series.  Figures  12a  and  12b  show  the  Time  Residual 
versus  T  for  a  transfer  angle  of  60  degrees.  Plots  for  other  transfer  angles  are  in 
Appendix  C.  The  maximum  absolute  value  of  the  Time  Residual  versus  transfer  angle 
for  5  and  15  terms  is  shown  in  Figures  13a  and  13b.  For  |T|  <  .25  15  terms  provided 
two  to  four  orders  of  magnitude  greater  accuracy  than  5  terms.  Recall  that  using  the 
true  coefficients  resulted  in  a  nearly  constant  four  orders  of  magnitude  difference. 
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Figure  12a.  Time  Residual  vs.  T  for  a  Transfer  Angle  of  60°  Using  5 

Terms  of  the  R  Series,  Method  A 
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R  SERIES  (METHOD  A).  15  TERMS 
THETA  =  60  DEGREES 
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Figure  12b.  Time  Residual  vs.  T  for  a  Transfer  Angle  of  60°  Using 

15  Terms  of  the  R  Series,  Method  A 
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Figure  13a.  Maximum  Absolute  Value  of  the  Time  Residual  vs.  Trans¬ 

fer  Angle  Using  the  R  Series,  Method  A,  for 
T  -  -  .25  to  .25 
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Figure  13b.  Maximum  Absolute  Value  of  the  Time  Residual  vs.  Trans¬ 

fer  Angle  Using  the  R  Series,  Method  A,  for 
T-  -.50  to  .50 
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Method  B:  R  Series  Coefficients  as  a  Function  of  Transfer  Angle 

Employing  the  least  squares  analysis  an  eighth  order  polynomial  was  found  to  best 
fit  the  plot  of  a,  versus  transfer  angle.  Therefore,  an  eighth  order  polynomial  was  then 
fit  to  the  plots  of  the  remaining  15  coefficients.  The  variance  increased  slightly  with 
each  successive  coefficient.  Figures  14a  -  14d  display  the  plots  for  the  first  four  coeffi¬ 
cients.  The  plots  for  the  fifth  coefficient  to  the  sixteenth  are  in  Appendix  D.  The 
Coefficient  Residual  (the  absolute  value  of  the  difference  between  the  true  and  the 
approximate  coefficient)  versus  transfer  angle  is  plotted  in  Figure  15.  Each  curve  has 
nearly  the  same  sinusoidal  shape.  The  distance  between  adjacent  curves  decreases  with 
each  successive  coefficient.  See  Appendix  D  for  a  plot  showing  Coefficient  Residual  for 
the  fifth  through  the  eighth  coefficients. 
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THETA  (RAD) 


Figure  14d.  a4vs.  Transfer  Angle 


The  accuracy  of  this  approximation  was  analyzed  as  before.  Figures  16a  and  16b 
show  the  Time  Residual  versus  T  for  a  transfer  angle  of  60  degrees.  Plots  for  other 
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transfer  angles  are  in  Appendix  E.  Figures  17a  and  17b  show  the  maximum  absolute 
value  of  the  Time  Residual  versus  transfer  angle  for  5  and  15  terms.  The  difference  in 
accuracy  between  5  and  15  terms  is  not  as  great  using  this  method  as  it  is  using  the  true 
coefficients  or  Method  A.  For  some  transfer  angles  5  terms  is  more  accurate  than  15. 
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Figure  16a.  Time  Residual  vs.  T  for  a  Transfer  Angle  of  60°  Using  5 

Terms  of  the  R  Series,  Method  B 
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Figure  16b.  Time  Residual  vs.  T  for  a  Transfer  Angle  of  60°  Using 

15  Terms  of  the  R  Series,  Method  B 
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Figure  17a.  Maximum  Absolute  Value  of  the  Time  Residual  vs.  Trans 

fer  Angle  Using  the  R  Series,  Method  B,  for 
T  --.25  to  .25 
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Figure  17b.  Maximum  Absolute  Value  of  the  Time  Residual  vs.  Trans¬ 

fer  Angle  Using  the  R  Series,  Method  B,  for 
7*  ■  -  .50  to  .50 

Continued  Fraction  Expansions 

In  the  continued  fraction  expansion  of  the  RI  series  the  Pk  and  Qk  terms  quickly 

approached  zero  as  k  increased.  The  twelfth  convergent  was  the  highest  that  could  be 
calculated  before  a  singularity  was  reached.  The  Pk  and  Qk  terms  for  the  R  series 
quickly  increased  as  k  increased.  The  highest  convergent  that  could  be  calculated  was 
the  eighteenth.  Therefore,  the  fifth  and  tenth  convergents  were  selected  to  compare 
with  the  series  using  5  and  15  terms  for  consistency. 

Figures  18a  and  18b  show  the  plot  of  the  Time  Residual  versus  T  for  a  transfer 
angle  of  60  degrees  using  the  continued  fraction  expansion  of  the  Rl  series.  Plots  for 
other  transfer  angles  are  in  Appendix  F.  The  accuracy  is  much  better  for  long  elliptical 
transfers  (positive  T)  than  the  RI  series  itself.  Figures  19a  -  19c  show  the  maximum 
absolute  value  of  the  Time  Residual  versus  transfer  angle  for  the  fifth  and  tenth  conver¬ 
gents.  Up  to  about  250  degrees  there  is  very  little  difference  between  the  accuracy  of 
the  fifth  and  tenth  convergents  for  1 71  <  .25.  However,  for  transfer  angles  greater  than 
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250  degrees  the  tenth  convergent  is  one  to  two  orders  of  magnitude  better.  For  \  t\  <  .50 
and  i T\  <  .85  the  tenth  convergent  is  more  accurate  than  the  fifth.  However,  for  some 
transfer  angles  the  difference  is  very  small  especially  over  the  interval  |7|  <  .85. 
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Figure  18a.  Time  Residual  vs.  T  for  a  Transfer  Angle  of  60°  Using 

the  5th  Convergent  of  the  C.  F.  Expansion  of  the  RI  Series 
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Figure  18b.  Time  Residual  vs.  T  for  a  Transfer  Angle  of  60°  Using 

the  10th  Convergent  of  the  C.  F.  Expansion  of  the  RI 
Series 
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Figure  19a.  Maximum  Absolute  Value  of  the  Time  Residual  vs.  Trans¬ 

fer  Angle  Using  the  C.  F.  Expansion  of  the  RI  Series  for 
T  m  -  .25  to  .25 
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Figure  19b.  Maximum  Absolute  Value  of  the  Time  Residual  vs.  Trans 

fer  Angle  Using  the  C.  F.  Expansion  of  the  RI  Series  for 
T  -  -  .50  to  .50 
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Figure  19c.  Maximum  Absolute  Value  of  the  Time  Residual  vs.  Trans¬ 

fer  Angle  Using  the  C.  F.  Expansion  of  the  RI  Series  for 
T  m  -  .85  to  .85 


Figures  20a  and  20b  show  the  plots  of  the  Time  Residual  versus  T  for  a  transfer 
angle  of  60  degrees  using  the  fifth  and  tenth  convergents  of  the  continued  fraction 
expansion  of  the  R  series.  Plots  for  other  transfer  angles  are  in  Appendix  G.  The  tenth 
convergent  provides  good  accuracy  for  long  elliptical  transfers  while  the  accuracy  of  the 
fifth  convergent  begins  to  fall  off  at  approximately  T  -  .6.  Figures  21a  -  21c  show  the 
maximum  absolute  value  of  the  Time  Residual  versus  transfer  angle  for  the  fifth  and 
tenth  convergents.  For  |T|  <  .25  and  |7|  <  .50  the  tenth  convergent  provides  one  to  two 
orders  of  magnitude  greater  accuracy  in  most  cases.  However,  for  transfer  angles  less 
than  or  equal  to  30  degrees  there  is  very  little  difference  over  the  interval  |T|  <  .25.  For 
1 7"  i  <  .85  the  accuracy  of  the  tenth  convergent  is  two  to  three  orders  of  magnitude  better 
for  transfer  angles  less  than  90  degrees.  At  about  230  degrees  the  fifth  convergent  is 
actually  more  accurate  than  the  tenth. 
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Figure  20a.  Time  Residual  vs.  T  for  a  Transfer  Angle  of  60°  Using 

the  ith  Convergent  of  the  C.  F.  Expansion  of  the  R  Series 
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Figure  20b.  Time  Residual  vs.  T  for  a  Transfer  Angle  of  60°  Using 

the  10th  Convergent  of  the  C.  F.  Expansion  of  the  R 
Series 
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Figure  21a.  Maximum  Absolute  Value  of  the  Time  Residual  vs.  Trans¬ 

fer  Angle  Using  the  C.  F.  Expansion  of  the  R  Series  for 
T  --.25  to  .25 
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Figure  21b.  Maximum  Absolute  Value  of  the  Time  Residual  vs.  Trans¬ 

fer  Angle  Using  the  C.  F.  Expansion  of  the  R  Series  for 
T  --.SO  to  .50 
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Figure  21c.  Maximum  Absolute  Value  of  the  Time  Residual  vs.  Trans¬ 

fer  Angle  Using  the  C.  F.  Expansion  of  the  R  Series  for 
7 --.85  to  .85 
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Comparison  of  Methods 

The  RI  series  is  more  accurate  than  the  R  series  by  two  to  three  orders  of  magni¬ 
tude  when  using  the  same  number  of  terms  as  shown  by  the  plots  in  Figures  22a  and 
22b.  The  R  series  requires  many  more  terms  be  used  in  order  to  obtain  the  same 
accuracy  as  the  RI  series. 
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Figure  22a.  Maximum  Absolute  Value  of  the  Time  Residual  vs.  Trans¬ 

fer  Angle  Using  5  Terms  of  the  R  Series  and  the  RI  Series 

for  T  -  -  .25  to  .25 


Using  the  alternate  methods  to  calculate  the  R  series  coefficients  did  not  introduce 
a  large  amount  of  error  when  5  terms  were  used.  The  accuracy  of  Method  A  was  nearly 
identical  to  the  accuracy  when  the  true  coefficients  were  used.  This  is  shown  in  Figure 
23a.  The  accuracy  of  the  alternate  methods  deteriorated  when  more  terms  of  the  series 
were  used  as  shown  in  Figure  23b. 
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Figure  22b.  Maximum  Absolute  Value  of  the  Time  Residual  vs.  Trans¬ 

fer  Angle  Using  15  Terms  of  the  R  Series  and  the  RI 
Series  for  T  -  -  .50  to  .50 


Figure  23a.  Maximum  Absolute  Value  of  the  Time  Residual  vs.  Trans¬ 

fer  Angle  Using  5  Terms  of  the  R  Series  With  the  True 
Coefficients,  Method  A,  and  Method  B  for 
T  -  -  .25  to  .25 
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Figure  23b.  Maximum  Absolute  Value  of  the  Time  Residual  vs.  Trans¬ 

fer  Angle  Using  15  Terms  of  the  R  Series  With  the  True 
Coefficients,  Method  A,  and  Method  B  for 
T  -  -  .50  to  .50 


Figure  24a  shows  that  the  fifth  convergent  of  the  continued  fraction  expansion  of 
the  RI  series  is  more  accurate  than  using  5  terms  of  e  series  for  transfer  angles  less 
than  260  degrees.  From  260  to  360  degrees  the  accuracy  of  each  is  approximately  equal. 
Figure  24b  illustrates  the  efficiency  of  the  continued  fraction  expansion.  The  tenth 
convergent  requires  eleven  coefficients  and  provides  nearly  the  same  accuracy  as  15 
terms  of  the  RI  series. 
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Figure  24a.  Maximum  Absolute  Value  of  the  Time  Residual  vs.  Trans 

fer  Angle  Using  5  Terms  of  the  RI  Series  and  the  5th 
Convergent  of  the  C.  F.  Expansion  of  the  RI  Series  for 
7  -  -  .85  to  .85 
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Figure  24b.  Maximum  Absolute  Value  of  the  Time  Residual  vs.  Trans 

fer  Angle  Using  15  Terms  of  the  RI  Series  and  the  10th 
Convergent  of  the  C.  F.  Expansion  of  the  RI  Series  for 
7- -.85  to  .85 


47 


The  continued  fraction  expansion  of  the  R  series  provided  approximately  three 
orders  of  magnitude  greater  accuracy  than  the  R  series  using  five  coefficients  for  each. 
This  is  shown  in  Figure  25a.  Figure  25b  indicates  the  continued  fraction  expansion  still 
provides  approximately  three  orders  of  magnitude  greater  accuracy  with  the  tenth  con¬ 
vergent  compared  to  15  terms  of  the  series. 
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Figure  25a.  Maximum  Absolute  Value  of  the  Time  Residual  vs.  Trans¬ 

fer  Angle  Using  5  Terms  of  the  R  Series  and  the  5th 
Convergent  of  the  C.  F.  Expansion  of  the  R  Series  for 
T  --.25  to  .25 


Using  the  continued  fraction  expansion,  the  R  series  approached  the  accuracy  of 
the  RI  series.  For  the  fifth  convergent,  the  accuracy  of  the  R  series  is  nearly  identical 
to  the  RI  series  for  transfer  angles  greater  than  250  degrees.  Below  250  degrees  the  RI 
series  is  still  approximately  an  order  of  magnitude  more  accurate.  This  is  shown  in 
Figures  26a  and  26b  for  two  time  intervals.  When  the  tenth  convergent  is  used,  the 
accuracy  of  the  two  series  is  nearly  identical  for  T  -  -  .85  to  .85  and  is  identical  for 
T  -  -  .25  to  .25  as  shown  in  Figures  27a  and  27b. 
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T  -  -.50  to  .50 
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Figure  25b.  Maximum  Absolute  Value  of  the  Time  Residual  vs.  Trans¬ 

fer  Angle  Using  15  Terms  of  the  R  Series  and  the  10th 
Convergent  of  the  C.  F.  Expansion  of  the  R  Series  for 
T  --.50  to  .50 


0  30  60  90  120  150  180  210  240  270  300  330  360 

ANGLE  (DEGREES) 


Figure  26a.  Maximum  Absolute  Value  of  the  Time  Residual  vs.  Trans¬ 

fer  Angle  Using  the  5th  Convergent  of  the  C.  F.  Expan¬ 
sions  of  the  R  and  the  RI  Series  for  T  -  -  .25  to  .25 
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0  30  60  90  120  150  180  210  240  270  300  330  360 

ANGLE  (DEGREES) 

Figure  26b.  Maximum  Absolute  Value  of  the  Time  Residual  vs.  Trans¬ 

fer  Angle  Using  the  5th  Convergent  of  the  C.  F.  Expan¬ 
sions  of  the  R  and  the  RI  Series  for  T  -  -  .85  to  .85 
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Figure  27a.  Maximum  Absolute  Value  of  the  Time  Residual  vs.  Trans¬ 

fer  Angle  Using  the  10th  Convergent  of  the  C.  F.  Expan¬ 
sions  of  the  R  and  the  RI  Series  for  T  -  -  .25  to  .25 
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MAX  ABS(TIME  RESIDUAL) 
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Figure  27b.  Maximum  Absolute  Value  of  the  Time  Residual  vs.  Trans¬ 

fer  Angle  Using  the  10th  Convergent  of  the  C.  F.  Expan¬ 
sions  of  the  R  and  the  RI  Series  for  T  «  -  .85  to  .85 
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IV.  Conclusions 


Performing  a  series  reversion  on  Eq  (5)  avoids  the  singularity  in  the  RI  series. 
However,  the  R  series  is  not  as  accurate  as  the  RI  series  with  the  same  number  of  terms. 
Another  disadvantage  of  the  R  series  involves  calculating  the  coefficients  which  is  very 
computer  intensive. 

Alternate  methods  of  calculating  the  R  series  coefficients  using  polynomial  approx¬ 
imations  were  found  to  be  nearly  as  accurate  as  the  true  R  series  coefficients.  However, 
these  methods  are  empirical  and  do  not  provide  any  advantage. 

Expanding  each  series  into  a  continued  fraction  provided  the  most  accurate  method 
of  calculating  the  semi-major  axis.  The  continued  fraction  expansions  were  most  not¬ 
ably  more  accurate  than  each  series  for  long  elliptical  transfers  approaching  the  mini¬ 
mum  energy  transfer  time.  The  R  series  continued  fraction  expansion  was  as  accurate  as 
the  RI  series  continued  fraction  expansion  with  the  benefit  of  no  singularities.  For  a 
specified  degree  of  accuracy  the  continued  fraction  expansion  is  more  advantageous  than 
the  series  because  fewer  series  coefficients  are  required. 

The  accuracy  of  each  of  the  methods  presented  was  dependent  upon  the  number  of 
series  terms  or  the  order  of  the  convergent,  the  transfer  time,  and  especially  the  transfer 
angle.  The  best  accuracy  obtained  in  this  study  was  a  Time  Residual  on  the  order  of 
10  ''.  For  a  given  number  of  terms  or  convergent  this  varied  by  as  much  as  several 
orders  of  magnitude  depending  upon  the  transfer  time  and  transfer  angle.  The  largest 
reduction  in  accuracy  usually  occurred  for  transfer  angles  between  270  and  350  degrees. 
As  T  approached  the  limit  for  which  the  number  of  terms  or  convergent  was  able  to 
provide  good  results  the  accuracy  was  very  sensitive  to  the  transfer  angle.  For  transfer 
times  well  within  the  limit  the  accuracy  did  not  vary  by  more  than  one  order  of  magni¬ 
tude. 
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V.  Recommendations  for  Further  Study 


During  the  course  of  this  research  the  radius  vectors  were  held  constant  while  the 
transfer  angle  was  varied.  Further  analysis  of  the  accuracy  may  be  accomplished  while 
varying  the  radius  vectors  and  the  transfer  angle.  Also,  further  analysis  of  the  accuracy 
may  reveal  the  optimum  number  of  terms  to  use. 

In  order  to  take  full  advantage  of  the  R  series,  a  more  efficient  method  of  obtain¬ 
ing  the  coefficients  must  be  found.  Further  investigation  of  the  R  matrix  may  reveal  a 
pattern  that  would  be  beneficial  in  order  to  generate  the  entire  matrix. 

While  the  continued  fraction  expansion  may  appear  to  be  an  attractive  method  of 
calculating  the  semi-major  axis  it  has  the  disadvantage  of  being  computer  intensive. 
Therefore,  it  may  be  desirable  to  approximate  the  continued  fraction  using  a  Chebyshev 
polynomial. 
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Appendix  A.  Time  Residual  vs.  T  Plots  for  Rl  Series 
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Figure  28a.  Time  Residual  vs.  T  for  a  Transfer  Angle  of  150 
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Figure  28b.  Time  Residual  vs.  T  for  a  Transfer  Angle  of  150 
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Appendix  B.  Time  Residual  vs.  T  Plots  for  R  Series 
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52a.  Time  Residual  vs.  T  for  a  Transfer  Angle  of  240 
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Figure  32b. 


Time  Residual  vs.  T  for  a  Transfer  Angle  of  240°  Using 
15  Terms  of  the  R  Series 
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53a.  Time  Residual  vs.  T  for  a  Transfer  Angle  of  330°  Using 

5  Terms  of  the  R  Series 
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Figure  33b. 


Time  Residual  vs.  T  for  a  Transfer  Angle  of  330°  Using 
i5  Terms  of  the  R  Series 


Appendix  C.  Time  Residual  vs.  T  Plots  for  Method  A 
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Figure  34a.  Time  Residual  vs.  T  for  a  Transfer  Angle  of  150°  Using 

5  Terms  of  the  R  Series,  Method  A 
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Figure  34b.  Time  Residual  vs.  T  for  a  Transfer  Angle  of  150°  Using 

15  Terms  of  the  R  Series,  Method  A 
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Figure  35a.  Time  Residual  vs.  T  for  a  Transfer  Angle  of  240°  Using 

5  Terms  of  the  R  Series,  Method  A 
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Figure  35b.  Time  Residual  vs.  T  for  a  Transfer  Angle  of  240°  Using 

15  Terms  of  the  R  Series,  Method  A 
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Figure  36a.  Time  Residual  vs.  T  for  a  Transfer  Angle  of  330“  Using 

5  Terms  of  the  R  Series,  Method  A 
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Figure  36b.  Time  Residual  vs.  T  for  a  Transfer  Angle  of  330°  Using 

15  Terms  of  the  R  Series,  Method  A 
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Figure  37d.  a„vs.  Transfer  Angle 
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Appendix  E.  Time  Residual  vs.  T  Plots  for  Method  B 
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Figure  39a.  Time  Residual  vs.  T  for  a  Transfer  Angle  of  l  SO"  Using 

5  Terms  of  the  R  Series,  Method  B 
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Figure  39b.  Time  Residual  vs.  T  for  a  Transfer  Angle  of  150°  Using 

15  Terms  of  the  R  Series,  Method  B 
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Time  Residual  vs.  T  for  a  Transfer  Angle  of  240°  Using 
5  Terms  of  the  R  Series,  Method  B 
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Figure  40b. 


Time  Residual  vs.  T  for  a  Transfer  Angle  of  240°  Using 
15  Terms  of  the  R  Series,  Method  B 
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lla.  Time  Residual  vs.  T  for  a  Transfer  Angle  of  330°  Using 

5  Terms  of  the  R  Series,  Method  B 


1.2E-001 
1.0E-001 
8.0E-002 
6.0E-002 
4.0E-002 
2.0E-002 
6.9E-01 8 
-2.0E-002 


R  SERIES  (METHOD  B).  15  TERMS 
THETA  =  330  DEGREES 


-4.0E— 002  1  i  '  '  i  |l  . . .  I  I  I  If  I"  I  "I1  I  |  I  I  I  'I  |  I  I  i  i  rrm  | 

-1.00  -0.75  -0.50  -0.25  0.00  0.25  0.50  0.75  1.00 

TIME  (t/tp  -  1) 


Figure  41b. 


Time  Residual  vs.  T  for  a  Transfer  Angle  of  330*  Using 
15  Terms  of  the  R  Series,  Method  B 


Appendix  F.  Time  Residual  vs.  T  Plots  for  C.  F.  Expansion  of  RI  Senes 
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Figure  42a.  Time  Residual  vs.  T  for  a  Transfer  Angle  of  150°  Using 

the  5th  Convergent  of  the  C.  F.  Expansion  of  the  RI  Series 
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Figure  42b.  Time  Residual  vs.  T  for  a  Transfer  Angle  of  150°  Using 

the  10th  Convergent  of  the  C.  F.  Expansion  of  the  RI 
Series 
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Figure  43&.  Time  Residual  vs.  T  for  a  Transfer  Angle  of  240°  Using 

the  5th  Convergent  of  the  C.  F.  Expansion  of  the  RI  Series 
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Figure  43b.  Time  Residual  vs.  T  for  a  Transfer  Angle  of  240°  Using 

the  10th  Convergent  of  the  C.  F.  Expansion  of  the  P»I 
Series 
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Figure  44a.  Time  Residual  vs.  T  for  a  Transfer  Angle  of  330°  Using 

the  5th  Convergent  of  the  C.  F.  Expansion  of  the  RI  Series 
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Figure  44b.  Time  Residual  vs.  T  for  a  Transfer  Angle  of  330°  Using 

the  10th  Convergent  of  the  C.  F.  Expansion  of  the  RI 
Series 
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Appendix  G.  Time  Residual  vs.  T  Plots  for  C.  F.  Expansion  of  R  Series 
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Figure  45a.  Time  Residual  vs.  T  for  a  Transfer  Angle  of  150°  Using 

the  5th  Convergent  of  the  C.  F.  Expansion  of  the  R  Series 
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Figure  45b.  Time  Residual  vs.  T  for  a  Transfer  Angle  of  150°  Using 

the  10th  Convergent  of  the  C.  F.  Expansion  of  the  R 
Series 
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Figure  46a.  Time  Residual  vs.  T  for  a  Transfer  Angle  of  240°  Using 

the  5th  Convergent  of  the  C.  F.  Expansion  of  the  R  Series 
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Figure  46b.  Time  Residual  vs.  T  for  a  Transfer  Angle  of  240°  Using 

the  10th  Convergent  of  the  C.  F.  Expansion  of  the  R 
Series 
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Figure  47a.  Time  Residual  vs.  T  for  a  Transfer  Angle  of  330°  Using 

the  5th  Convergent  of  the  C.  F.  Expansion  of  the  R  Series 
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Figure  47b.  Time  Residual  vs.  T  for  a  Transfer  Angle  of  330°  Using 

the  10th  Convergent  of  the  C.  F.  Expansion  of  the  R 
Series 
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